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Abstract 

The free energy density of the two-dimensional 0(3) non-linear cr-model is cal- 
culated at finite temperature and finite spatial extent. We make both an analytic 
calculation in the perturbative regime and a Monte-Carlo study at low tempera- 
tures. We show that using the fixed-point action instead of the standard Wilson 
action leads to a great reduction of the cut-off effects. 
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1 Introduction 



The calculation of thermodynamical quantities in Yang-Mills gauge theories 
using lattice regularization suffers from huge cut-off effects. Decreasing the lat- 
tice spacing a, the cost of such calculations grows with - , where in SU(3) 
gauge theory e is around fO. This requires a very careful analysis to remove 
the cut-off effects. Therefore in field theory one looks for improved actions 
which give good results even on lattices with bad resolution. Wilson's idea of 
the renormalization group (RG) [1,2] aims at that direction : one integrates 
out the short-distance degrees of freedom and replaces them by a smaller 
set of effective degrees of freedom, which describe long-distance physics. One 
possibility is to use a blocking procedure introduced originally by Kadanoff. 
Theoretically there exists a perfect lattice action, that predictions are abso- 
lutely free of any cut-off effects. We use here the fixed-point (FP) action, 
which is the classical perfect action and is expected to be a good approxima- 
tion for the perfect action. For the 0(3) non-linear cr-model, the model we are 
considering here, the FP action was developed in [3]. Recently in a nice work 
[4], Alessandro Papa pointed out the great advantage of the FP action for 
SU(3) gauge theory: The free energy density near the phase transition point 
showed an impressive reduction of cut-off effects compared to earlier results 
derived by other lattice actions [5] . However the few points calculated in that 
paper do not allow a complete view of the behaviour in the decrease of cut- 
off effects. In the case of the a-model we can continue this project and give 
the free energy density over a large temperature region. Hence we check the 
cut-off effects arising from different lattice actions both in the high tempera- 
ture regime (perturbation theory) and at low temperatures using Monte-Carlo 
(MC) simulations. 

The Euclidean continuum version of the partition function for the 0(3) non 
linear a-model in two dimensions reads as follows: 



Jds(x) 



e -M(S) 



where S is a three component vector with S 2 (x) = 1 and measure 

DS(x) = l[d 3 S(x)d(S 2 (x)-l) . (2) 



The term /3A(S) is the continuum action with periodic boundary conditions 
for the field S at finite temperature T = l/L t and spatial extent L s . 

PA cont (S) = I jdxojdx, d,S(x)d,S(x) (3) 
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The free energy density f is given by the relation 

Z = e~ fV 



(4) 



where V is L t x L s and f depends on both L t and L s . To get the temperature 
and spatial size dependent, physical part of this quantity we subtract the 
contribution resulting from the square box, in which we set L t equal to L s . In 
the perturbative regime and for infinite spatial extent we get — f/T 2 = n/3 
which is the two-dimensional case of an ideal gas of massless scalar particles. 

This letter consists of the following parts: in section 2 we give a short review 
of the ideas which lead to the FP action in the cr-model. The next section 
presents the perturbative calculation of the free energy density and a compar- 
ison between different lattice actions. Section 4 treats the cluster-MC analysis 
necessary for the calculation in the non-perturbative, low temperature regime 
and shows the corresponding results. 



2 FP action 

We summarize now some important results about the FP action of the 0(3) 
a-model which we use in this letter. For a detailed discussion we refer the 
reader to ref. [3]. 

We divide our original lattice into 2x2 blocks labeled by n B - The four spins 
S n inside a block form, by taking some mean, a new spin variable R nB on a 
lattice with a lattice unit twice as large as the original one. The new action 
A!{R) is defined by integration over the original lattice: 



where T is the kernel of the RG transformation with the normalization that 
ensures the invariance of the partition function 





(6) 



In the limit (3 — > oo the blocking kernel T has the form 



T(R,S) = kJ2( E Sn ~Rn B - E S: 



n 



) 



(7) 
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where k (fixed here to 2) can be tuned to find a short-ranged FP-action. In 
the limit mentioned above equation (5) reduces to a saddle point problem, 
giving 



A'(R) = mm [a(S) + T(R, S) 



{S} 

and for the FP action 

A FP {R) = minj^pp (S) + T(R, S)\ . (9) 
{s} I J 

For parametrizing it, one can write in general 
Afp(S) = -\^p(j-){l - S n S n+r ) 

n,r 

+ c(n 1 ,n 2 ,n 3 ,n 4 )(l- S n3 S n4 ) + ...,(10) 

ni,ri2,n3,n4 

where it is a significant help, that the first two functions p and c can be 
calculated analytically. The Fourier transform of p has the form 



P(q) ±L (q + 2vr/) 2 to (ft/2 + ^) 2 3« 



The summation is over the integer vector I = (l , / x ) and (q + 27r/) 2 = (q + 
2nl ) 2 + (qi + 27r/!) 2 . Furthermore one can show that p fulfills the FP equation 
for a free scalar theory with a gaussian block-transformation. 

Beside the two functions p and c the resulting couplings were found in [3] 
by a numerical fitting procedure and recently improved in [6] to control the 
topological effects. We use the parametrization obtained in [6]. It consists of 
24 different couplings which can be put on a 2 x 2 lattice. When working with 
such an action to measure the free energy density one expects to see cut-off 
effects, due to the following three points: 

• The FP action is not the perfect action. 

• In the simulation a parametrized form of the FP action is used which intro- 
duces some parametrization errors. 

• Finally, on extremely small lattices (N t = 2) the action is disturbed by the 
small size, which causes cut-off effects that decrease exponentially with N t . 
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3 Perturbative results 



In a normal perturbative treatment of eq. (1), one introduces a Faddeev-Popov 
term to avoid the zero mode problem [7]. One arrives finally, when changing 
to Fourier space, at the following form of the partition function: 



where the sum here just leaves away the unwanted zero mode k = (0,0), 
Ylk = Z)fc,fc^(o,o)- A s a remainder of the Faddeev-Popov term we have to keep 
the logarithm of the volume in the partition function. Due to equation (3) the 
Fourier sum above is discreet and one can use any regularization scheme to 
extract the continuum value f con t- When using lattice regularization one must 
replace the continuum propagator k 2 in eq. (12) with the corresponding lattice 
propagator, for example with the function p of eq. (11) in the FP case. We 
then check the cut-off effects by considering the ratios fw /fcont, fs /fcont, fa/ fcont 
and fpp/fcont, where the abbreviations "W, S " denote the use of Wilson and 
Symanzik action, whereas "fi" means an action with a nearest neighbour and 
a diagonal coupling, whose ratio is the same as in the FP action. Although 
there is no problem to give f as a function of £, £ = L • T — L s / L t , we present 
the cut-off effects for ( = 3. Table 1 gives the results f CO nt/fiatt for different 
lattice schemes, a different number of points N t in the time direction (i.e. for 
different resolutions) and N s = 3 x N t . An approximation of the decrease of 
the cut-off effects can then be given through a function of the form: 

1 + '^ + "'W + c '^ + --- ■ (13) 



The Wilson action shows the expected quadratic decrease (fig. 2a), whereas 
in the Symanzik case (fig. 2b) the fit is consistent with a = 0. Figure 2c 
shows the behaviour of f n and figure 2d finally presents a perfect exponential 
decrease of the cut-off effects for the FP action. The smooth curve is given by 
1 + A ■ exp(—BN t ), where A and B are constants. 



4 MC methods 



In a MC simulation the partition function cannot be measured directly, on the 
other hand, using equations (1,4), the derivative of the logarithm of Z with 
respect to (3 gives: 

^ In Z = -Z- 1 [DS A(S) e^ A ^ . (14) 



5 



Table 1 

Numerical values for different lattice actions with N s = 3 x N t 

# N t fw/fcont fcont/fs fcont/ffi fcont/fFP 



2 


1.1606 


0.964058 


1.016701 


1.006314 


4 


1.0455 


1.001507 


1.004622 


1.000016 


6 


1.0182 


1.000155 


1.001773 


1.000000 


8 


1.0098 


1.000024 


1.000729 


1.000000 


10 


1.0062 


1.000007 


1.000240 


1.000000 



where (At) denotes the expectation value of the action on the lattice N t x 
N s and (*4. vac ) the same on the square lattice which we subtract as vacuum 
contribution. 

In our simulations we used a cluster algorithm [8], which almost completely 
eliminates the critical slowing-down, which was especially important for the 
measurement of the mass gap. For the measurement of the action we performed 
8 x 10 5 sweeps per /3-value and estimated the error by the method of bunching. 
The size of the error is of the order of the circles in figure 2, where one 
finds the signal of the integrand in eq. (15). The smooth curve is a spline fit 
interpolation, which we used for the integration, where f3 had to be chosen so 
small that the signal for that value disappeares. To estimate the error of this 
integration we compared the value resulting from the spline curve with the 
value one gets, joining the single points in fig. 2 by straight lines. This error is 
not bigger than half the size of the symbols in figures 3, 4 and 5. To compare 
the cut-off effects we measured the action for Wilson and FP parametrization. 

The above procedure gives just the free energy density as a function of (3. To 
extract physical results we must relate every /3-value to a temperature. This 
can be done by measuring another physical quantity, the mass gap. In a lattice 
simulation one can only get the dimensionless quantity m-a, where the lattice 
spacing a is a function of (3. Using the relation T = l/N t a we find: 



One obtains 




(15) 



T 



1 



(16) 



m maN t 
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Fig. 1. Figures a-d show the disappearence of the cut-off effects for different lattice 
actions. The smooth curves are fitted functions explained in the text. 



where N t is the number of points in the time direction of the lattice used for 
the measurement of (.At)- We then extracted the quantity m ■ a(f3) for some 
/3-values by taking a square lattice which was about seven times larger than 
the correlation length to avoid finite size effects. Inserting these values in eq. 
(16) we got a finite number of points T/m {(3). To our data we added further 
results of the mass gap in the FP case from ref. [6], while for the Wilson action 
we used the results from [9]. 

We present first (fig. 3) a comparison between Wilson action and FP action. 
Fig. 4 gives the detailed results obtained with the FP action and the last figure 
(fig. 5) shows a spline interpolation of the FP data. 
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Fig. 2. The signal of the integrand in eq. (15) for a 4 x 12 lattice is plotted against the 
values of (3. The smooth curve is the spline interpolation we used for the integration. 



5 Summary 

The most important result of this work is the great reduction of cut-off effects 
with the FP action compared to other lattice actions. In the perturbative 
regime the cut-off effects disappear exponentially fast and are nearly lost on 
a lattice as small as N t = 4. Concerning the cut-off effects, our results of MC 
simulations show within the range of statistical error a similar behaviour in 
the low-temperature region. Although we are working with a relatively small 
statistical error we are not able to detect any power-law in the decrease of the 
cut-off effects. 
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Fig. 3. The discreet points give the free energy density — f/T 2 for different lattice 
resolutions with Wilson action, whereas the smooth curve is a spline interpolation 
of the data resulting from a 3 x 9 lattice with FP action. The error of the free 
energy density is smaller than the size of the symbols. The horizontal line shows 
the perturbative value for ( = 3 in the high temperature limit. 
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Fig. 5. The smooth curves show a spline interpolation of the data for 2 x 6, 3 x 9 
and 4 x 12 lattices in the FP case (fig. 4). Within statistical error the values for the 
5 x 15 lattice lie on the curve of the 4 x 12 lattice. 



